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Abstract 

Chiral perturbation theory gives direct and unambiguous predictions for 
the form of various two-point hadronic correlators at low momentum in terms 
of a finite set of couplings in a chiral Lagrangian. In this paper we study 
the feasibility of extracting the couplings in the chiral Lagrangian (through 
1-loop order) by fitting two-point correlators computed in lattice QCD to the 
predicted chiral form. The correlators are computed using a pseudofermion 
technique yielding all-point quark propagators which allows the computation 
of the full four-momentum transform of the two-point functions to be obtained 
without sacrificing any of the physical content of the unquenched gauge config¬ 
urations used. Results are given for an ensemble of dynamical configurations 
generated using the truncated determinant algorithm on a large coarse lattice. 
We also present a new analysis of finite volume effects based on a finite volume 
dimensional regularization scheme which preserves the power-counting rules 
appropriate for a chiral Lagrangian. 
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1 Introduction 


The introduction of the concept of an effective chiral Lagrangian |l|, and the subse¬ 
quent development of systematic chiral perturbation theory [Q, have led to remark¬ 
able clarihcations in our understanding of the properties of quantum chromodynamics 
(QCD) in the low-momentum regime. Much of the phenomenology of the interac¬ 
tions of hadrons at low energies can be understood and organized using chiral theory, 
where the basic degrees of freedom are colorless light hadrons, rather than the quarks 
and gluons appropriate for a description at (presumably arbitrarily) short distance 
scales. Of course, the parameters in a chiral effective Lagrangian are ultimately 
determined by the underlying microscopic theory, i.e. the QCD Lagrangian. The 
connection between the underlying theory and the chiral Lagrangian is however com¬ 
pletely nonperturbative, so any determination of chiral parameters from QCD must 
have recourse to a systematic nonperturbative calculational procedure: at the present 
time, this means lattice gauge theory. 

Perhaps the most direct and unambiguous predictions of chiral perturbation the¬ 
ory refer to the low-momentum behavior of various two-point correlators of hadronic 
densities and currents. Typically, the low-momentum expansion of these correlators 
is ordered in a series of increasing powers of the momentum-squared with the 
power-counting rule rri^ ~ p^. The leading order term is then determined by the 
lowest-dimension part £2 of the full chiral Lagrangian £chir = £2 + £4 + ••• treated 
at tree level, the next to leading term at low-momentum in hadronic correlators is 
obtained from use of £2 at one-loop level, together with £4 at tree-level, and so on. 
In this paper we show that the accurate extraction of up to the leading three terms of 
this expansion is perfectly feasible by simulation of the relevant correlators in lattice 
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gauge theory. We use exclusively unquenched conhgurations throughout, so there are 
no issues of quenched chiral logs, for example. However, the great expense needed 
to generate decorrelated dynamical conhgurations at light quark mass (our conhgu¬ 
rations were generated using the truncated determinant algorithm (TDA) [§) makes 
it obviously desirable to extract the maximum possible physical content from each 
gauge conhguration. Consequently, we have used a stochastic pseudofermion tech¬ 
nique to obtain all-point quark propagators for each gauge conhguration. This 

means that the Fourier transform used to determine each two-point correlator studied 

A(p) = ^ < 0{x)0{y) > (1) 

contains a double sum over all points on the lattice and therefore a factor of the 
lattice volume more terms than the equivalent correlator studied with single-source 
or single-sink quark propagators. 

Another aspect of calculations performed on discrete space-time lattices is the 
necessary presence of an infrared (large-distance) cutoh. In this paper we present 
results for hadronic correlators on a physically large but coarse lattice (6“^ lattice 
with lattice spacing a=0.4 F, O(a^) improved gauge action) at two diherent sea- 
quark masses (note: in all the calculations the valence and sea-quark masses are kept 
the same). Using a hnite volume generalization of dimensional regularization, we 
have computed the appropriate hnite volume modihcations of the chiral formulas of 
Ref.p). Of course, the systematic errors induced by lattice discretization will require 
availabity of a variety of unquenched lattices at dihering couplings. 

The presentation of our results in this paper is organized as follows. In Section 
2 , we review the all-point pseudofermion technique i,§ used to obtain quark prop¬ 
agators for the measurement of hadronic correlators (for a more detailed study of 
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the statistical features and computational cost of this method see |0). In Section 
3, we present results for the pseudoscalar-pseudoscalar two-point function for our 
unquenched lattice ensemble. In particular, we show that consistent hts to the pre¬ 
dicted chiral form can be obtained, determining the corresponding one-loop chiral 
Lagrangian parameters to a (statistical) accuracy of a few percent. Moreover, the 
htted higher order coefficients are quite small, suggesting that chiral perturbation 
theory can continue to be accurate for surprisingly high momentum (the chiral hts 
are extended up to ~2.5 GeV^). In fact, the calculations presented here of higher 
order chiral coefficients really provide a quantitative nonperturbative basis for esti¬ 
mating the accuracy of low order chiral perturbation theory. In Section 4, we repeat 
this procedure for the axial vector-axial vector current correlator. Section 5 contains 
a discussion of the hnite volume dimensional regularization technique we have used to 
derive systematically the hnite volume corrections to the chiral formulas for hadronic 
correlators. Of course other important systematic errors, primarily those due to lat¬ 
tice discretization, can only be addressed once a much larger selection of dynamical 
lattice ensembles at various lattice spacings and with improved quark actions are 
available. 
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2 Computing Hadronic Correlators with All-point 


Quark Propagators 

The dynamical configurations used in the present work represent a very high invest¬ 
ment in computational effort. Each new configuration generated using the truncated 
determinant algorithm requires evaluation of several hundred low eigenvalues of 
the hermitian Wilson-Dirac (or clover) operator. On a 1.5 GHz Pentium 4 proces¬ 
sor, this takes 3 minutes for a 6'^ lattice and about 80 minutes for a 10^x20 lattice. 
Decorrelation of physical quantities can take anywhere from 50 to several hundred 
such updates. The very high cost of generating properly decorrelated unquenched 
gauge configurations makes it essential to squeeze the maximum physical information 
content from each available configuration. Conventional quark propagators obtained 
by (say) a conjugate gradient algorithm only provide the quark propagation ampli¬ 
tude from a single source vector (either point or smeared) to any other point on the 
lattice: accordingly, the computation of a two-point correlator such as Eq(l) would 
require fixing either point x or y (typically, at the lattice origin), which obviously 
sacrifices a great deal of physical information in the gauge configuration. The results 
presented in this paper will involve correlators computed from all-point quark prop¬ 
agators which give the quark propagation amplitude from any point on the lattice to 
any other point. 

Here we describe briefly an approach (originally suggested by Michael and Peisa in 
the context of static quark systems, f^) to obtaining such propagators by simulating 
bosonic pseudofermion fields. For each quark propagator needed in the hadronic 
observable of interest (in the case of two-point current correlators, there are two 
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quark propagators which must be multiplied and traced appropriately) introduce a 
bosonic pseudofermion held (p^a with action (m a lattice site, a the spin-color index, 
Q the Wilson or clover operator): 


S{(j)) = (p^Q'^Qcj) 

= H = ^,Q = H^ 


( 2 ) 

(3) 


For a given hxed background gauge held A, simulating the pseudofermion held with 
the preceding action produces the following correlator (<< O » means the average 
of O relative to the measure e~^) : 


(pma^Pnb 


(4) 


(-/? )ma,nb 

(5) 


(,Q T5)ma,nfe 

(6) 


Note that separate pseudofermion helds are needed for each quark propagator as av¬ 
erages of four-point bosonic pseudofermion amplitudes will produce contractions with 
the wrong sign relative to the corresponding fermionic 4 quark amplitudes. However, 
it turns out that the computational ehort required is quite manageable, allowing us to 
obtain sufficiently accnrate all-point propagators with only a few times the compnter 
time required for a conventional conjngate-gradient evaluation of the corresponding 
single-sonrce propagator. The reasons for this are twofold: 

(i) The pseudofermion average << .... >> is efficiently implemented by a heat-bath 
npdate of pseudofermion helds. As the psendofermion action is a Gaussian one, the 
statistical properties of this simulation are essentially trivial (gauge conhguration up¬ 
dates, by contrast, involve an action with a highly nonlinear dependence on the held 
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degrees of freedom), and contain no surprises. 

(ii) For a fixed gauge field, most quantities decorrelate after a few pseudofermion 
sweeps. For low-momentum quantities, the presence of low eigenmodes of the hermi- 
tian Wilson operator H can lead to long autocorrelations, but these can be handled 
by projecting out the corresponding low modes, which reduces the condition number 
and restores a more rapid decorrelation (see [^). However, even without such projec¬ 
tions, we are able to obtain the momentum dependence of the correlators sufficiently 
accurately to allow good fits to the chiral behavior. 

It is easy to see that the computation of multipoint hadronic correlators involv¬ 
ing n quark propagators can be reduced to convolutions of n pseudofermion fields, 
rapidly computed by fast Fourier transform (FFT). (This is so both for local and 
smeared hadronic operators, although the quantities studied in this paper are exclu¬ 
sively local densities and currents.) For example, the full 4-momentum transform 

of fho 2-point pseudoscalar correlator is given by 

A(x,?/) = < 0|T{T(a;)75T(a;) T(|/)75T(2/)}|0 > 

= - < tr((Q“^75),,y(Q“^75)yx) > 

'y H)xa 

ab 

= - « » (7) 

which becomes an easily evaluated fast Fourier transform of products of pseud¬ 
ofermion fields: 

A(g2) = - « FFT(xtH0)(q)FFT(0tHx)(-q) » (8) 

Our main purpose in this paper is to show that the chiral behavior of QCD 
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hadronic correlators can be studied directly in momentum space and higher order 
chiral couplings extracted with high statistical accuracy. The unquenched lattices 
used here lattices are a large ensemble (800 conhgurations) of physically large, coarse 
lattices recently used in a study of stringbreaking 0 (lattice spacing 0.4 Fermi, but 
with O(a^) improvement of the gauge action 0: thus, the gauge action includes 
the usual plaquette term as well as a twisted rectangle (“trt”) operator), and at a 
sea-quark kappa value of 0.2050. Quark eigenmodes up to 420 MeV are included 
exactly in the determinant. In this lattice ensemble, the determinant contribution 
to the measure corresponds to two degenerate light quark flavors (corresponding to 

~200 MeV, where the lattice scale is hxed from stringtension measurements). 

We conclude this section by describing the computational effort required for eval¬ 
uation of the pseudofermion average (7) on a single gauge conhguration. For the 6^ 
lattices, a single heat-bath update of the two pseudofermion helds 0, y requires 0.366 
sec. on a 1.5 GHz Pentium-4 processor. The convolutions and FFT operations re¬ 
quired to obtain the desired four-momentum held A(q'^) in (8) require an additional 
0.024 sec. and are performed after every 2 heat-bath updates of 0, y. The hnal pseud¬ 
ofermion average for A{q‘^) was obtained from 20000 measurements, corresponding to 


2.1 Pentium-4 hrs. 


3 Extracting Chiral Parameters from Pseu¬ 
doscalar Correlators 

Among the most basic predictions of chiral perturbation theory are expressions for 
the low momentum behavior of correlators of hadronic densities and currents, hrst 
derived by Gasser and Leutwyler @]. In this section we study the extent to which 
a lattice computation of the two-point function of the isovector pseudoscalar density 
V’A 75 t/’ (here, Tj are the isospin generators for SU(2), ip the isodoublet quark field) 
can be used to constrain the parameters of the chiral effective Lagrangian. We work 
throughout with two flavors of degenerate light dynamical quarks, so the relevant 
chiral group is SU(2)xSU(2), and we will adhere as much as possible to the notation 
of Ref 0, the results of which are in any event restricted to exactly this two-flavor 
case. 

Dehning t± = |(ti ± ir 2 ), the correlator computed in (7-8) corresponds to 

Aps-ps{q^) = J d^x < V5r_75?/>(a;)Vir+75V’(l/) > (9) 

with the low-momentum chiral behavior (in infinite volume Euclidean space) 0: 

Aps-PsU/) ^ - fti) + 0(«t) (10) 

and with the pseudoscalar decay constant defined as 

G^6ij =< Ol'tjjTi'y^iplTTj > ( 11 ) 

The quark condensate < 0\'pi'pj\0 >= F^B6ij hxes the constant B, and G, hi are 
couplings in the O(p^) chiral Lagrangian. Our object here is to determine the latter 
as accurately as possible using lattice data. The fits will also provide information on 
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the size of even higher order terms (e.g. O(p^) in the psendoscalar correlator), which 
is clearly relevant to the issne of accnracy of chiral perturbation theory through one 
loop order. 

We begin by describing the results obtained from our ensemble of 800 6^ lat¬ 
tices. These large coarse lattice conhgurations are O(a^) improved with respect to 
the gauge action, but the quark action is unimproved Wilson, so the results obtained 
will necessarily contain large systematic effects (relative to the corresponding con¬ 
tinuum parameters), although we shall see that the statistical errors are extremely 
small. The conhgurations were generated using the TDA algorithm P] to include 
virtual quark-loop effects of a doublet of degenerate light sea quarks exactly up to a 
quark off-shellness of about 420 MeV. The average value obtained for Aps- Ps{q^) for 
this ensemble is shown in Fig.l. For nonzero momentum, the statistical errors are 
smaller than the symbol size. The range of corresponds to lattice values from 0 to 
16, with a unit of q^ on the lattice corresponding to 0.25 GeV^ in physical units. The 
comparatively large error at the zero momentum point is attributable to the autocor¬ 
relations induced by low eigenmodes of H (see 0 for a solution of this problem). In 
fact, our results for the nonzero momentum modes are sufficiently accurate that we 
will discard the zero-momentum point entirely in performing the hts to Eq(lO). 

As the psendoscalar correlator (9) involves by dehnition only local operators, it 
may be wondered whether there is much hope of extracting an accurate pion mass 
in (10). The excited state contamination which would normally require the use 
of smeared sources/sinks to obtain an accurate pion mass appears as contributions to 
the higher order chiral contributions (of order 0(g°, q^...)). A conventional cosh £t 

to 800 smeared-local correlators at k =0.2050 gives an alternative determination of 
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Table 1: Dependence of fitted pion mass on momentum fitting range. 


g^range 


xVd-o.f 

1-8 

0.324± 0.033 

3.6 

1-9 

0.352T 0.026 

2.7 

1-10 

0.422T 0.020 

2.6 

1-11 

0.449± 0.017 

2.3 

1-12 

0.480T 0.015 

2.2 

1-13 

0.530± 0.013 

2.7 


the pion mass for this ensemble: = 0.396 ± 0.007 (or about 200 MeV in physical 

units, with the scale determined from string tension measurements [0: the use of an 
unimproved Wilson action on this coarse a lattice makes spin-dependent scales such 
as the rho unreliable). The inevitable breakdown of chiral perturbation theory in 
the ultraviolet means that we cannot expect the £t to give meaningful results if the 
momentum range is too large. Indeed, we hnd that, allowing the pion mass to vary 
freely in a £t of the form 

^ 2 I 2 + ^2 + (12) 

the £t values for increase well beyond the physical value if the fit range is extended 
beyond q'^ ~2.5 GeV^ (see Table 1). There is a broad minimum in the chisquared per 
degree of freedom for the momentum range 1 < < 10 — 12: choosing the UV cutoff 

in our fits at = 10 (lattice units,~ 2.5Gev^) gives a fitted pion mass very close to 
the value found using smeared source correlators. 

A five parameter £t of the formula (12) performed on the range 1 < g^ < 10 yields 
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Figure 2: Fit of measured pseudoscalar correlator A(g^), 1 < < 10 

(see Fig. 2) the following results: 

= 1.52 ±0.029 
= 0.422 ±0.020 
A 2 = 0.622 ± 0.010 
A 3 = -0.0460 ±0.0015 

A 4 = 0.00163 ±0.00007 (13) 

with a x^/d.o.f = 2.6. The pion mass is consistent within errors with the value 
quoted above, determined from a study of smeared-local correlators. The remarkable 
statistical accuracy of this method is apparent in the errors for the leading order 
residue Ai and the one-loop chiral parameter A 2 , both of which are determined at 
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the 2% level for the statistical error! Moreover, the higher order coefficients q‘^) 

terms) are very small, which bodes well for the accuracy of chiral perturbation theory 
through one-loop, even at quite high momenta. Repeating the £t with the pion mass 
held hxed at the value M^- =0.396 determined from hts of smeared-local correlators 
gives a somewhat improved x^/d.o.f=2.2, with the htted parameters statistically 
indistinguishable from the results in (13). However, removing the term in (12) 
results in hts with much higher y^/d.o.f. >5, and £t parameters Ai,A 2 and 
differing substantially from the values obtained in (13) if the momentum range is 
extended beyond q^ =6, corresponding to 1.5 GeV^ in physical units. The value of 
the pion decay constant will be obtained from a study of the axial-vector correlator in 
the next section: we find TV =0.187±0.011. The critical kappa value is most readily 
extracted from the predicted topological charge distributions of Leutwyler and Smilga 


lb| : one hnds fi:c=0.2067, corresponding to a bare quark mass for n = 0.2050 of 


m =0.020. On the other hand, using the value for Ft, quoted above, the bare quark 
mass can be extracted from the chiral Ward identity 


m = 


FttM^ 


(14) 


where from (10,13) we find G-,, =1.74±0.02. Accordingly (taking Mj, =0.422), (14) 
yields a bare quark mass 

m = 0.0191 ±0.0021 (15) 


consistent with the Kc value extracted from the topological analysis. 
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4 Extracting Chiral Parameters from Axial Cur¬ 


rent Correlators 


In this section we shall describe our results for the nonperturbative lattice evalua¬ 
tion of the axial vector correlator, using the same lattice ensemble described in the 
preceding section. Defining 

^Ax-Axi^i) = Jd^x< (le) 

we have used the pseudofermion technique to extract the contracted scalar quantity 
(in Euclidean space) 


^AX-Ax{q ) = x9ij,u^'aX-Ax{(1) 


(17) 


The prediction of chiral perturbation theory for this contracted axial-vector correlator 
can be summarized in the formula 10 : 


A.«-.4x(9^) - 4F; + - h,)i- + 0(,‘) 


( 18 ) 


where the term involving — h 2 represents a contribution from the O(p^) chiral 
Lagrangian. The constant term —in (18) arises from a contact term oc (5^(a; — y) 
in the coordinate space correlator which is highly UV-divergent and not accessible 
from our lattice calculation, so we have allowed this constant to float freely in our 
fits. 

We have computed A ax-ax { 9 “^) using the pseudofermion approach: for the en¬ 
semble of 800 lattice conhgurations, the results are shown in Fig. 3. Again, the 
statistical errors are quite small. With the htting formula 


AAx-Axiq ) = 


2 r^A/f2 + 


(19) 
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Figure 3: Measured axial-vector correlator ^ax-ax{.(f)) fit 0 < < 10 

a conservative fit using only the range from 0 to 4 (lattice units) and setting the 
higher order term to zero, one finds =0.213±0.010 with a y^/dof of 0.8/2. 
In this £t we have fixed the pion mass at the value 0.422 determined in the preced¬ 
ing section from the pseudoscalar correlator hts. The higher order coefficients are 
determined as Ai =0.8157±0.0040, A 2 =0.0602±0.0004 respectively. The statistical 
error on these higher order chiral coefficients is again very small, indicating the high 
statistical content available in an all-point approach. The contact term Ai is stable 
to three signihcant figures with respect to changing the htting range to 0 < <3 

or 0 < g^ <5, as it is completely dominated by ultraviolet contributions (see discus¬ 
sion above). Another point worth noting here is that our momentum-space approach 
allows the extraction of meaningful chiral parameters at light quark mass from a local- 
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local axial-vector correlator, something which would be impossible using conventional 
Euclidean coordinate space hts. 

Keeping only terms up to order in the chiral htting formula, one hnds in the case 
of the axial-axial correlator that the chisquared deteriorates much more rapidly as 
higher momenta are included in the htting range than in the case of the pseudoscalar 
correlators discussed in the preceding section. If a term is included in the ht, then 
the htting range can be extended considerably. In Fig. 3 the chiral ht obtained for 
a htting range 0 < g^ <10 is shown: the y^/dof for this ht is 14/7. The extracted 
chiral parameters in this case are: 

= 0.187 ±0.011 
Ai = 0.8148 ±0.0038 
A 2 = 0.0777 ±0.0003 

ylg = -0.00442 ± 0.00002 (20) 

5 Finite Volume Effects 

5.1 Dimensional Regularization at Finite Volume 

In order to reliably extract the parameters of the chiral Lagrangian out of equations 
(12) and (19) and the values in equation (13) and (20), we need a generalization of 
equations (10) and (18) that includes hnite volume ehects. At the same time this 
generalization should preserve the power counting rules appropriate for chiral pertur¬ 
bation theory. Here we present such a generalization by extending the dimensional 
regularization scheme to hnite volume with periodic boundary conditions. We assume 
our system lives in a box with the momentum taking values where is 
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an integer. 


We start by converting massive propagators to exponential form using Schwinger 
parameters: 

, 2 ^ 2 = ( 21 ) 

Jo 

and replacing the inhnite volume loop integrals over k (where k is a. d-dimensional 
vector) by hnite sums 

d^k 


A-d 


/i 




( 22 ) 


Once the k-dependence is exponentiated as in (21), we will need to evaluate sums 


like 


E e-"(—) =0FJ(—A) 


nu=—oo 


c? 


so that (22) will involve where X is the function 

1 - 1-00 

1(A) ^ ^ i: 


TT 


(23) 


(24) 


We will also need to evaluate a generalization of X(A): 

1 - 1-00 


TT 


(25) 


where we shall need both X and J in the physically interesting limit A —a 0 corre¬ 
sponding to large but hnite volume. 

Using the Poisson sum formula one can easily write equations (24) and (25) in a 
form useful to take that limit: 


and 


1 H-CX) 2 2 1 2 

_ , , L X-> TT m L , TT 

1 {\) = — Y. e A ~^(l + 2e A + 

V A ra=—oQ V A 


1 2 2 
L p_A ^^ _ TT m / /n \ 

J{\,p) = —/=^ E ^ ^ cos(7rpm) 

V A 


m=—oo 


(26) 


(27) 
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These equations are all we need to develop the generalization of dimensional regular¬ 
ization to finite volume that preserves the power counting rules appropriate for chiral 
perturbation theory. 

Consider for example the quadratically divergent integral 

•-■'■‘-iSfwh’ 

where we work at finite volume. Using the procedure outlined above and changing 
variables we get 

^2 = r V(A) (29) 

4 Jo 

Replacing the expression (26) for X we can easily check that the first term in the 
expansion is identically equal to the infinite volume integral while the other terms 
correspond to ultraviolet finite exponentially small finite volume corrections. The 
divergence as d —4 is therefore equal to the infinite volume divergence and the MS 
or MS prescriptions remain unmodified. 

If mC in the exponent above is large, the next to leading term in the expansion 
of X^ approximates well the finite volume correction of 5*2. We can make d = 4 since 
the term is finite, getting 

2 2 
_ m C _TT_ 

^next to leading 2 ^ ( r\ (oi\\ 

where Ki [x) is a modified Bessel function. This is exponentially small when mC gets 
large. Further corrections fall exponentially at an even faster rate. For the lattice 
ensemble studied previously, m£~2.5 so this approximation is inadequate and (29) 
must be evaluated more carefully (see below). 

As another example of the finite volume dimensional regularization procedure 
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consider the relation 


for the quartically divergent integral. For inhnite volume this is the only consistent 
dehnition, as there is no available dimensionful quantity other than /i, which can 
only appear logarithmically in the hnite volume part at d = 4. At hnite volume this 
argument does not work because we have now a new scale £. However, as we show 
next, the relation is also valid at hnite volume. 

Writing 1 in the “integrand” as {m? + and introducing again the 

Schwinger parameter we can easily transform 5*4 into 

1 . 4 —d.jj-d/2 nQo O 

= d\e-\x-^)n\) (32) 

where x = as usual. Separating X‘^(A) into and X'^(A) — A'^/^ and replacing 
each part into the above equation, we get for the hrst part 

..4—rf_c//2 «oo f) 

^divergent ^ ^ _ ^)A^/2 (33) 

JO o\ 

which is identical to the inhnite volume version of 5*4 and therefore it must be zero 
by the standard argument. In fact, analytically continuing equation (33) to d < 0 
and integrating by parts one can see that it is zero as it must be. 

The non trivial result is that the hnite volume correction 

1.4 d^dj‘1 /.Qo 

^convergent ^ ^ _ ^d/2) (34) 

Lx Jo UA 

is also explicitly zero for any d, as one can see by replacing the expansion (26) for 
X(A) above and integrating by parts. 

It is now clear that with the procedure outlined in equations (21) to (27) we 
can compute the hnite volume ehects in the dimensional regularization scheme. In 
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the next subsection we shall employ this technology to compute the hnite volume 
corrections to the two-point hadronic correlators. This will allow us to estimate 
quantitatively the systematic errors in the results of Sections 3,4 due to hnite size. 


5.2 Finite Volume Corrections to Two-Point Hadronic Cor¬ 
relators 

Let us hnd hnite volume corrections to the pseudoscalar correlator and the axial 
vector correlators using the dimensional techniques described above. First, we note 
that the hnite-volume dependence of correlators calculated through O(p^) in chiral 
perturbation theory arises from the one-loop integrals using the lowest order chiral 


Lagrangian for the chiral vertices (see Ref |^). Since one-loop graphs come from £ 2 , 


the evaluation of pseudoscalar and axial vector correlators involves loop integrals of 
only the type (28) @]. So in calculating the hnite volume corrections to and 

Fjr to one loop order we can simply add the hnite correction to this loop integral to 
the mass logarithm terms in those constants obtained in Ref [@. 

The Gtt; and modihcations are then given by 


and 


where 


G. 

_^ ^inf. vol. F qP 


(36) 

M^- 

► (AC'™') + 


(36) 

Fn 

T-iinf. vol. ^ 


(37) 

2 “ 4 £ 2 „ 

/ dXe J^ A - 

'0 


(38) 
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which corresponds to the hnite volume part of S 2 of equation (29). Since the coeffi¬ 
cient in the exponent in the integrand, M^£^/47r^, is small(~ 0.16), the contribution 
from large A is substantial. So the original form (24) for X(A) is used. The numer¬ 
ical integration using Mathematica yields 5^=0.002200. To hnd the hnite volume 
corrections to the leading order, the values of the constants F7r=0.187±0.011, and 
M.,r=0.422±0.020 found in the previous sections are used for F, and M. With 
above, the hnite volume correction to is 0.007 and XA -0.012. In the hnite volume 
term of G^there is a constant B (practically, the quark condensate) which was not 
extracted directly from the previous hts. To the leading order = 2BF ||^. Thus 
the hnite volume correction to Gj^ becomes —Gt^S^ /2X^ to this order. Using a central 
value of G7r=1.74, we hnd -0.05 for the hnite volume correction to Gjr. We see that 
in all cases the hnite volume corrections are on the order of a few percent even for 
the fairly light quark mass (Mjr ~ 200 MeV) used in the dynamical lattices studied 
in Sections 3,4. Of course, these lattices were physically large (2.4 F*^): on smaller 
lattices for light quark masses, the hnite volume corrections will be more important. 
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